Energy Decay in Josephson Qubits from Non-equilibrium Quasiparticles 
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We calculate the energy decay rate of Josephson qubits and superconducting resonators from 
non-equilibrium quasiparticles. The decay rates from experiments are shown to be consistent with 
predictions based on a prior measurement of the quasiparticle density n qp = 10 //an 3 , which sug- 
gests that non-equilibrium quasiparticles are an important decoherence source for Josephson qubits. 
Calculations of the energy-decay and diffusion of quasiparticles also indicate that prior engineered 
gap and trap structures, which reduce the density of quasiparticles, should be redesigned to improve 
their efficacy. This model also explains a striking feature in Josephson qubits and resonators - a 
small reduction in decay rate with increasing temperature. 

PACS numbers: 03.65.Yz, 74.50+r, 74.25.Nf 



Superconducting integrated circuits are a leading can- 
didate for scalable quantum information processing [lj], 
with recent experiments showing accurate state con- 
trol of coupled circuits @. Josephson qubits are rela- 
tively straightforward to couple because their relatively 
large size enables quantum interactions via simple wiring. 
However, these same wires also make the qubit suscep- 
tible to various decoherence mechanisms, and thus a de- 
tailed understanding of all loss mechanisms is critically 
needed 0]. Non-equilibrium quasiparticles are known 
to be an important decoherence mechanism for charge 
qubits, since quasiparticles hopping on and off the qubit 
island produce an unpredictable change in state 0, [||. 
How do quasiparticles affect the coherence of charge- 
insensitive qubits, such as phased, flux, and transmonQ 
devices? 

In this letter, we present a simple and general model 
that describes how non-equilibrium quasiparticles pro- 
duce energy decay in superconducting qubits and res- 
onators. We relate predictions to the quasiparticle den- 
sity n qp : based on its prior measurement [8[ , we calculate 
an energy decay rate Ti that could possibly explain cur- 
rent limits. Although other decay mechanisms, such as 
surface dielectric loss[9j, are important, experimentalists 
need a detailed theory of all mechanisms in order to prop- 
erly design experiments that separate out and measure 
energy decay, and to ultimately improve qubit perfor- 
mance. We conclude that quasiparticles must be consid- 
ered in present experiments. 

Our model is described in four parts. First, we calcu- 
late the energy decay rate for a Josephson qubit, using 
a generalization of the environmental P(E) theorvfioj]. 
Since the prediction is based on tunneling of non- 
equilibrium quasiparticles, we next calculate their den- 
sity versus energy from standard electron-phonon rates, 
showing the injection energy is unimportant. Third, we 
estimate quasiparticle diffusion lengths and comment on 
the design of effective quasiparticle traps. Finally, we 
consider and discuss possible sources of non-equilibrium 



quasiparticles, using energy estimates. 

No equilibrium quasiparticles are expected in a super- 
conductor well below the critical temperature T <C T c , 
because of the exponential suppression of excitations 
from the superconducting gap A. However, an exper- 
iment on Cooper-pair boxes has measured a significant 
quasiparticle density n qp ~ 10 //im 3 , arising from an un- 
known source Q. Once excited, quasiparticles can readily 
tunnel through a Josephson junction. 

Charge associated with quasiparticles allows a tun- 
neling event to couple energy between the quasiparticle 
and the qubit, as described by the environmental P(E) 
theory [10| - When a qubit state of energy E w decays, its 
energy is added to the final quasiparticle state. Energy 
is transferred, however, in only a fraction of the tunnel- 
ing events for the case of large junction capacitance C 
that we consider here. For the simple case of the transfer 
of charge q, the probability for qubit decay from a single 
tunneling event is p ~ (q 2 /2C)/E W . We extend the P(E) 
theory to account for quasiparticle tunneling having both 
electron-like (q — — e) and hole-like (q — e) c harg e trans- 



fer, which gives a qubit decay probability [111 [1 
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where u and v are the (energy dependent) BCS amplitude 
factors for the left (I) and right (r) junction electrodes. 

For a tunnel junction with normal state resistance Rt, 
the total qubit decay rate is obtained by summing over 
all possible quasiparticles and accounting for filled and 
empty states. For the tunneling of quasiparticles from 
left to right, we integrate over all possible initial energies 
E to obtain a rate 
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where pi, r (E) — E/(E 2 — A 2 ) 1 / 2 is the normalized quasi- 
particle density of states, and // jT . is the non-equilibrium 
occupation probability. The factors u r , v r , p r and f r 
are computed for the final quasiparticle state with en- 
ergy E + Eio, and the coherence factor is opposite in 
sign as for charging effects^. The total decay rate Ti 
is obtained by also summing the quasiparticle tunneling 
rate from right to left T i . We note the similarity of 
this formula to the Mattis-Bardeen theory for metallic 
conductors 13L 14 1. 



The total density of non-equilibrium quasiparticles is 
given by n qp = 2D{E F ) J p{E)f(E) dE, where D(E F ) is 
the density of electron states at the Fermi energy. The 
integral for the decay rate can be solved for the physically 
relevant case of contributions from f(E) only near the 



gap, as motivated below, 
ratetla is 



For small n qp , the total decay 
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The decay rate scales as the fractional quasiparticle oc- 
cupation, which is normalized by the total density of 
Cooper pairs D(E F )A. 

We can immediately check the relevance of this pre- 
diction by plugging in typical parameters for qubits. 



In Ref. [8j, a quasiparticle density of 10 //zm 3 was ex- 
perimentally measured in charge qubits. For the alu- 
minum superconductor that is typically used, we have 
A/k = 2.1 K and D(E F )A = 2.8 • 10 6 /^m 3 . Typical 
phase qubit parameters are E\o/h — 6 GHz, R F — 200 fi, 
and C = 1 pF, which yields from Eq. ([4]) a qubit decay 
rate T\ = 1/(2.1 /is). This rate is reasonably close to 
what is measured for our qubit [16] T% ~ 550ns, espe- 
cially considering Ref. @ possibly had lower n qp because 
of quasiparticle traps. 

Non-equilibrium quasiparticles also produce energy de- 
cay in resonators from loss in superconducting wires. Fol- 
lowing the calculations in Ref. 17j , the resonator quality 
factor Q can be expressed as a function of the quasi- 
particle density Here, the bulk conductivity of the 
superconductor is given by the Mattis-Bardeen theory, 
from which the surface impedance can be expressed in 
various limits. Geometric factors expressing the non- 
uniform current density in the coplanar transmission line 
are also included, to yield in the thick-film limit 
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Here, A ~ (iTpoO'nA/h)^ 1 / 2 is the penetration depth cor- 
responding to the surface kinetic inductance po\, a n is 
the normal state conductivity, s is the width of the cen- 
ter conductor, g = 1.0 and g m = 0.31 are geometrical 
factors expressing the effect of non-uniform current den- 
sity on the kinetic and magnetic inductance, respectively, 



and 7 = 1/2 is a factor concerning the interface. Us- 
ing s — 3 /im, A — 50 nm for Al, a resonator frequency 
of 6 GHz, and quasiparticle parameters described previ- 
ously, we compute Q ~ 10 7 , which corresponds to a decay 
time of 300 ps. This time is a significantly longer than 
found above for qubits, and agrees with experimental ob- 
servation that decay times are longer for resonators [3] . 

A more detailed understanding of non-equilibrium 
quasiparticles can be obtained by calculating their occu- 
pation probability f(E) using electron-phonon scattering 
rates[lj|- We first consider a bulk superconductor. The 
lifetime of a quasiparticle at energy e to scatter to any 
energy e' by emitting or absorbing a phonon of energy 
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where To « 400 ns is the measured characteristic electron- 
phonon time for Al. The lifetime of a quasiparticle state 
with energy e to recombine with another quasiparticle of 
any energy e' is 
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Note the similarity in this formula to Eq. ([6]) , but here the 
energy of the emitted phonon is greater than 2 A. With 
the term /(e'), this rate is proportional to the density 
of quasiparticles, implying that the recombination rate 
slows down once the density of quasiparticles becomes 
small. 

Equations ([6]) and ([7]) are used to solve for the steady- 
state occupation probability f{E) of quasiparticles for 
a given (non-equilibrium) injection rate r qp . With the 
recombination rate being proportional to the square of 
the density of quasiparticles, we solve this nonlinear dif- 
ferential equation numerically. Because the quasiparticle 
relaxation rate at high energy is dominated by electron- 
phonon scattering, which is number conserving for quasi- 
particles, the exact energies at which the quasiparti- 
cles are injected are unimportant for the computed low- 
energy distribution. For simplicity, we inject quasiparti- 
cles with energies between 2.8 A and 3 A. 

We plot the quasiparticle occupation probability versus 
energy in Fig. [1] for an injection rate r qp /2D(E F )A = 
1.7 • 10~ 10 /t that gives a total quasiparticle density 
n qp /2D(E F )A — 1.8 • 10~ 6 , approximately the same as 
for Ref. Q. At phonon temperatures T > 170 mK, a 
thermal occupation fx of quasiparticles dominates. At 
T ~ 140 mK, an exponentially decaying occupation is 
calculated, but at an elevated occupation due to the 
non-equilibrium quasiparticles. For temperatures below 
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FIG. 1: Plot of quasiparticle occupation versus energy for 3 
phonon temperatures T p — 140, 70, and mK (top to bot- 
tom). An injection rate was chosen to match the experimental 
quasiparticle density 10 //im 3 , equivalent to n qp /2D(EF)A = 
/ fpdE = 1.8 • 10 -6 . The state occupation at energies shown 
does not depend on the injection energies. The gray line shows 
the thermal occupation fr at 140 mK. 



is D = 60w qp cm 2 /s, where v qp = (1 - A 2 /E 2 ) 1 / 2 is 
the normalized velocity, which depends on energy and 
approaches zero at the gap. Quasiparticles diffuse a dis- 
tance (Dt) 1 / 2 in time r; we plot this length in Fig. [3] 
for both the electron-phonon and electron-electron de- 
cay times. We estimate that the characteristic distance 
for quasiparticles to equilibrate are of millimeter lengths, 
much larger than the size of qubit devices. In normal 
metal traps, quasiparticles lose energy at much short 
length scales, of order 30 /im. Finally, quasiparticle en- 
ergy is ultimately removed from the superconductor by 
recombination and emission of a phonon with energy 
> 2 A. Because phonons are long-lived at low temper- 
ature, especially for crystalline substrates, they may bal- 
listically travel across the chip and have their energy 
redeposited as quasiparticles anywhere in the supercon- 
ductor. All of these estimates imply that a conservative 
design should have quasiparticle trapping structures of 
millimeter or greater size, not just local structures placed 
around the tunnel junctions. 

Although the density of quasiparticles depends on ex- 
perimental details of chip design, it is possible to make 
rough power estimates to help understand the generation 



about 70 mK, the curves become independent of tem- 
perature. This calculation is consistent with experi- 
ments on Cooper-pair-box experiments [20| . where non- 
equilibrium quasiparticle tunneling rates became inde- 
pendent of refrigerator temperature below about 70 mK. 

The total quasiparticle density at T = is numerically 
found to be n qp /2D{E F )A = 0.U{T O r qp /2D{E F )A) 1 / 2 . 
We find the computed decay rate is well approximated 
by Eq. g]). 

The calculated temperature dependence of the decay 
rate is plotted in Fig.[^a), which shows negligible change 
at low temperatures T < 40 mK, and then a small re- 
duction in rate up to T p = 120 mK. This unusual pre- 
diction of an improvement in decay rate with increasing 
temperature results from non-equilibrium quasiparticles 
having, on average, larger energy and thus smaller contri- 
butions from density of states and coherence factors. At 
higher temperatures, the thermal generation of quasipar- 
ticles rapidly increases the decay rate. This temperature 
dependence well explains data taken on superconduct- 
ing resonators [2lJ- Measurements of Ti(T) for a phase 
qubit [161] is shown in Fig.[^b), where a reduction and 
then an increase in decay rate is in qualitative agreement 
with predictions 22]. The simplicity of our model, which 
does not account for the small occupation of the qubit 
ground state, diffusion of quasiparticles, or gap inhomo- 
geneities, may explain the differences. 

Previous experiments have used engineered gap[23j| 
and trap 0, [24[ structures to reduce the density of quasi- 
particles. Are they large enough to be completely effec- 
tive? The quasiparticle diffusion constant for aluminum 
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FIG. 2: (a) Plot of calculated qubit energy relaxation rate 
versus phonon temperature T p for a quasiparticle density of 
n qp /2D(E F )A = 1.8 ■ 10~ 6 and E 10 /h = 6 GHz. This shape 
of the curve is relatively insensitive to changes in Eio and n qp . 
Note the small decrease in decay rate with increasing temper- 
ature up to 120 mK, and then a rapid increase due to ther- 
mally generated quasiparticles. (b) Experimental measure- 
ment of relaxation rate Ti for a phase qubit. The relaxation 
decreases, and then increases with temperature, in qualitative 
agreement with theory. The slow change at T p > 140 mK is 
probably due to non- idealities, such as the small occupation 
of the ground state (greater than 90% for T v < 100 mK, but 
only 77% at T p = 130 mK and 55% at T p = 160 mK). 



4 



"F io 1 






-i— i 


\ e-e 




D) 






§10° 




- 


: usion 




e-ph 


1"' 







1.1 1.2 1.3 1.4 

Quasiparticle energy E/A 



1.5 



FIG. 3: Plot of average quasiparticle diffusion length versus 
energy for both the electron-phonon f /V s (e-ph) and electron 
recombination 1/T r (e-e) decay times. The e-e interaction 
dominates for E < 1.1 A. Parameters are D = 60w qp cm 2 /s 
(Al), T p = 0, and n qp = W/fim 3 . 



mechanism. The model presented here found a quasipar- 
ticle injection rate r qp ~ 2.4 • 10 3 /s/um 3 . Assuming a 
total superconductor volume of 0.1 ^m x 10 mm 2 and an 
energy A per quasiparticle, the total power to the chip 
is 0.06 pW. 

The power load from various sources can be compared 
with this number. Cosmic rays have a flux of ~ 0.6/cm 2 s 
and deposit an energy of ~ lMeV/mm, yielding for a 
50 mm 2 chip of thickness 0.5 mm a power 0.02 pW. Back- 
ground radioactivity is typically of order that coming 
from cosmic rays. If coaxial lines allows ~ 4K thermal 
radiation to be transmitted with a bandwidth 100 GHz, 
a power ~ 5 pW is found. Blackbody radiation from 1 K 
gives about 60 pW of power load to the outside of the chip 
mount. Materials also slowly release energy at low tem- 
peratures, with amorphous Si02, polycrystalline Al, and 
teflon showing heat release of ~ lOOpW/gr, ~ 20pW/gr, 
and ~ 2pW/gr respectively that slowly decays with 
time 251. From these estimates, 



non-equilibrium quasi- 
particles are clearly plausible. 

In conclusion, we have developed a model that predicts 
qubit energy decay from the density of non-equilibrium 
quasiparticles. Predictions of the model are in reasonable 
agreement based on a prior experimental measurement 
and, at a minimum, show that this decoherence mech- 
anism should be carefully considered. The quasiparti- 
cle energy distribution calculated from simple electron- 
phonon interaction qualitatively agrees with measure- 
ments, and we argue that engineered gap and trap struc- 
tures should be of large (millimeter) size. With a good 
understanding of the basic physics, we believe experi- 



ments can be performed to more carefully test for non- 
equilibrium quasiparticles, and to ultimately improve the 
coherence of superconducting qubits. 
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